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Long-range quantum lattice systems often exhibit drastically different behavior than their short- 
range counterparts. In particular, because they do not satisfy the conditions for the Lieb-Robinson 
theorem, they need not have an emergent relativistic structure in the form of a light cone. Adopting 
a field-theoretic approach, we study the one-dimensional transverse-field Ising model with long-range 
interactions, and a fermionic model with long-range hopping and pairing terms, explore their critical 
and near-critical behavior, and characterize their response to local perturbations. We deduce the 
dynamic critical exponent, up to the two-loop order within the renormalization group theory, which 
we then use to characterize the emergent causal behavior. We show that beyond a critical value of 
the power-law exponent of the long-range couplings, the dynamics effectively becomes relativistic. 
Various other critical exponents describing correlations in the ground state, as well as deviations 
from a linear causal cone, are deduced for a wide range of the power-law exponent. 


I. Introduction and Summary 

Long-range interactions arise in a wide range of phys¬ 
ical systems. Examples include NV centers and other 
solid-state defects [1-4], excitons (Frenkel excitations) in 
organic solids [5], polarons [6], Shiba chains [7, 8], and 
photon-mediated interactions between superconducting 
qubits [9]. Furthermore, long-range interactions emerge 
very naturally—and are often unavoidable—in atomic, 
molecular, and optical (AMO) systems, for example van 
der Waals (1/i?®) interactions between Rydberg atoms 
[10, 11] or polaritons [12], magnetic or electric dipole- 
dipole (1/i?^) interactions between atoms or molecules 
[1-4, 10, 13-15], and variable-range (l/i?“) interactions 
between atoms in multimode cavities [16] or trapped ions 
[17-20]. While the quantum-critical behavior of short- 
range interacting models has been extensively studied, 
quantum criticality and phase transitions and their uni¬ 
versal properties are less fully explored in the presence of 
long-range interactions. 

An important difference between short- and long-range 
interacting systems, namely the nature of any emergent 
causal structure, can be characterized by their response 
to a local perturbation. The celebrated Lieb-Robinson 
bound demonstrates that even nonrelativistic quantum 
systems exhibit a linear ‘light cone’ bounding a causal 
region, outside of which the response function is expo¬ 
nentially suppressed [21], provided the interactions are 
short-ranged. This bound enforces the emergence of a 
‘relativistic’ causal behavior even in condensed matter 
systems. For long-range power-law interactions, on the 
other hand, the light cone may be sublinear, and the 
bounds on the influence of a local perturbation are much 
less stringent. In an early work by Hastings and Koma, 
the boundary of this region is shown to be at least log¬ 
arithmic rather than linear [22]. Recent works have fur¬ 
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ther explored the causal consequences of long-range in¬ 
teractions [23-27], and Ref. [28] has recently improved 
the Hastings-Koma bound by constraining the causal 
region algebraically rather than logarithmically. How¬ 
ever, tighter bounds are not ruled out, and it remains 
an open question whether a linear light cone emerges for 
generic power-law interacting models beyond some criti¬ 
cal power-law exponent. 

In this paper, we study the causal structure of one¬ 
dimensional long-range lattice models in the vicinity of 
a quantum critical point, and adopt a field-theoretical 
approach, complementary to the extensive literature on 
Lieb-Robinson-type bounds for long-range lattice mod¬ 
els. Specifically, we compute correlation functions and 
causal response functions in the ground states of contin¬ 
uum field theories governing the near-critical behavior 
of the long-range interacting transverse-field Ising model 
(TFIM), and also a fermionic model with long-range hop¬ 
ping and pairing terms. We identify the so-called dy¬ 
namic critical exponent that defines the relative scaling 
of space and time coordinates, which characterizes the 
causal structure of the underlying lattice model close to 
its critical point. We show that linear light cones emerge 
above critical values of the power-law exponent of the 
long-range couplings, which depend on the spin/fermion 
model as well as whether the model is at or away from 
criticality. In both cases, we also identify the critical 
exponent characterizing the decay of correlations in the 
ground state. Furthermore, it is shown in detail that the 
response to a local perturbation obtains a general scaling 
form, which is determined by the value of these two expo¬ 
nents. Note that approximate numerical approaches to 
many-body models with long-range interactions [29-32] 
also exist (in the context of ultracold systems or else¬ 
where). In contrast, our field-theory treatment is well 
suited to extracting universal aspects of the long-distance 
and long-time behavior of the many-body system, and 
specifically the critical exponents (beyond their mean- 
field values) in the thermodynamic limit, which are usu¬ 
ally difficult to access via numerical methods. 
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Table I: The exponents describing two long-range lattice models at criticality: TFIM and Fermion Model. The long-range 
interaction between Ising spins or long-range hopping and pairing terms in the fermionic model is assumed to fall off as 
~ where a > 0. The exponent 2 denotes the dynamic critical exponent; z = 1 defines a linear light cone, while z < 1 

corresponds to a sublinear causal region. The exponent 9 characterizes the decay of the correlation function in the ground state 
at criticality {9 being the anomalous dimension rj reported in Sec. Ill for the TFIM). For the TFIM, and in an intermediate 
range of the power-law exponent, the critical exponents should be computed in a series of epsilon-expansion with e = 3 (t/ 2 — 1 
and c((t) a rather complicated expression approximately given by (;{a) « 1/[24(1 + a‘^)]. Even away from the critical point, the 
correlation functions decay as a power law ~ (not shown in this table). 
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Table II: The causal behavior of long-range TFIM and Fermion Model at or away from criticality. For long-range models with 
a sufficiently rapidly decaying power-law, the causal behavior is described by a linear light cone. For exponents smaller than 
a critical value, which depends on the specific model, the causal behavior is not linear. This critical value of the power-law 
exponent, for the models at criticality, is given by the onset of the dynamic critical exponent deviating from 2 : = 1, cf. Tab. I. 
At criticality, the response functions Da (TFIM) and (Fermion Model) are described by the critical exponents a and 9, 
and the general scaling functions ga (TFIM) and Ho- (Fermion Model) as detailed in the text. For the critical TFIM model 
with a > 7/4, and the fermionic model with a > 1, the response function describes a linear light cone. Away from criticality, 
linear/non-linear behaviors are predicted by a simple analysis of the dispersion relation. Analytical expressions describing 
various regimes are provided in the text. For the noncritical Fermion Model, long-range hopping or pairing cases exhibit 
different causal behavior. 


It is important to note that Lieb-Robinson-type 
bounds are generally state-independent and agnostic to 
many details of the underlying model, requiring only a 
lattice with a finite-dimensional local Hilbert space on 
each site [21], and an interaction with some prescribed 
spatial decay. While these bounds are very general, 
naive applications to long-range interacting models re¬ 
sult in unphysically large causal regions, and obtaining 
even qualitatively tight rigorous bounds remains an open 
problem. 

In this paper, by limiting ourselves to ground states 
of specific models and local quenches on such ground 
states, we find a rich structure not captured in the more 
general Lieb-Robinson statements. Crucially, our results 
provide explicit examples of models with particular sub- 
linear causal regions. Thus, while existing Lieb-Robinson 


bounds constrain the speed with which information can 
propagate from above, our results place a lower limit 
on how far (towards lower speeds) such upper bounds 
can ultimately be pushed. One might expect that the 
availability of low-energy excitations at a quantum crit¬ 
ical point should facilitate faster-than-light propagation 
of the response to local perturbations in the presence of 
long-range couplings; we confirm this intuition by com¬ 
paring critical and noncritical regimes. In this sense, our 
study of critical points also lends some heuristic support 
to the conjecture that our results may actually coincide 
with the best possible Lieb-Robinson bounds for long- 
range interacting systems. However, this is certainly 
speculative; our treatment is complimentary to the Lieb- 
Robinson type approach and provides some intuition for 
how the most refined versions of such bounds might look. 
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but does not constitute a rigorous result on how such 
bounds must look. 

We also stress that our definition of a causal region 
and linear/non-linear causal behavior is not necessarily 
precisely equivalent to that used in the context of Lieb- 
Robinson bounds and information theory. In the latter 
context, the shape of the causal region is determined by 
requiring that the Lieb-Robinson bound falls below some 
threshold value outside of it, which naturally defines a 
region in space and time within which appreciably large 
signals can, in principle, be sent. In our definition, the 
light-cone shape is tied to the dynamic critical exponent, 
which is a natural and standard identification from the 
point of view of field theory. We show that this expo¬ 
nent controls how fast information reaches a given point 
by identifying it with the space-time scaling of the first 
local maximum of the response function, and thus there 
is clearly a close connection to the information-theoretic 
definition. We find this definition necessary to extract a 
light-cone shape from a specific model, since in the ab¬ 
sence of perfectly ballistic transport the height of this 
maximum generically decays in space and time, and con¬ 
stant contours do not extend asymptotically to large dis¬ 
tances and times. This issue presents itself for short- 
range interacting models as well, where in general the 
proper light cone shape is only obtained by ignoring the 
decay of the peak (local maximum) [33] , with the under¬ 
standing that the absence of perfectly ballistic dynamics 
is a model-specific phenomenon and should not influence 
the Lieb-Robinson bound. However, whether this notion 
exactly coincides with that of information theory is an 
interesting question in need of further investigation. 

We briefly discuss the methods and quantities of inter¬ 
est relevant to our investigation of critical and dynam¬ 
ical aspects of long-range interacting lattice models in 
Sec. II. We present our results for the long-range inter¬ 
acting TFIM in Sec. Ill, and for the fermionic model with 
long-range hopping and pairing terms in Sec. IV. We have 
studied in detail (i) the correlation/response functions, 
(ii) for spin/fermion models, (hi) at/away from critical¬ 
ity, and (iv) for different ranges of the exponent charac¬ 
terizing the power-law couplings. In all these models, we 
are always sufficiently close to criticality that the correla¬ 
tion length is large compared to the lattice spacing, and 
a continuum description is valid. For the benefit of the 
reader, we have summarized our main results in Tables I 
and II. We have performed a detailed renormalization- 
group (RG) calculation up to the two-loop order for the 
TFIM, and argued, on the basis of RG, that the criti¬ 
cal exponents obtained in the fermionic model are exact 
(i.e., mean-field exponents are exact, and would not re¬ 
ceive corrections from interactions). For both models, 
short-range interactions give rise to a dynamical expo¬ 
nent 0 = 1, which indicates a linear light cone and rel¬ 
ativistic dynamics. We explore in detail how sufficiently 
slow-decaying power-law couplings can give rise to sub- 
linear light cones, with 0 < 1. 


II. Methods and quantities of interest 

In this work, we rely heavily on scaling and renormal¬ 
ization group theory, which provides a systematic way to 
integrate out short-wavelength degrees of freedom (at the 
scale of lattice spacing, for example) in order to find an 
effective description of the physics at long wavelengths. 
We often find that general quantities of interest take a 
simple scaling form involving universal exponents that 
determine how fast correlations fall with distance, or how 
space and time coordinates scale with respect to each 
other, see Tables I and II. Such exponents can be approx¬ 
imately determined (at the mean field level) via a simple 
power counting, which is the first step of a systematic RG 
treatment. The knowledge of mean-field exponents can 
be then used to determine whether interactions affect the 
universal behavior of the system. In fact, it can be ar¬ 
gued, as we will often do, that many types of interaction 
become less and less important at long wavelengths (the 
corresponding coefficients are suppressed along the RG 
flow). In this case, they do not affect the universal behav¬ 
ior, and we say that interactions are not relevant in the 
sense of RG. In other cases where interactions are rele¬ 
vant, we systematically use renormalization group theory 
to determine the critical exponents beyond their mean- 
field values. This approach is particularly appealing since 
even a complicated theory can be efficiently described by 
a small set of exponents. 

In the remainder of this section we define the various 
quantities of interest reported in Tables I and II, separat¬ 
ing them into two subsections based on relevance to the 
TFIM or the fermionic model. We will focus on prop¬ 
erties of these quantities specifically in the limit of long 
times and large distances. 

TFIM. —In the model studied in Sec. Ill, we shall fo¬ 
cus on the 0 -component of the spin, Sf, that defines the 
Ising order parameter (i denoting the lattice site in a ID 
chain). We are primarily interested in universal prop¬ 
erties of correlation functions and causal response func¬ 
tions, 

iGAt,t-j) = {rs^{t)s;m, ( 1 ) 

iD,{t,t-j) = Q{t){[S!{t),S;{0)]), (2) 

where the operator S'f(t) is defined in the Heisenberg 
picture. The subscript a denotes the exponent of the 
long-range interaction potential V{i — j) ~ \/\i — 
with nearest-neighbor terms only, we have cr = 00 . In 
the above equations, 0 is the Heaviside step function, T 
is the time-ordering operator, and the expectation values 
are computed in the ground state of the system Hamilto¬ 
nian (see Sec. HI for specific details). We have also used 
the translation symmetry to write the two-point func¬ 
tions as a function of the distance (in units of lattice 
constant) between the two points and their time differ¬ 
ence. The correlation function probes the inherent corre¬ 
lations of the ground state, and is nonzero even at equal 
times t = 0. The response function characterizes causal- 
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ity in the system, namely how fast information propa¬ 
gates from a given point to another, and is constrained 
by the Lieb-Robinson bound to decay exponentially out¬ 
side of a linear light cone for short-range interactions [ 21 ], 
Fermion model .—In the fermionic model studied in 
Sec. IV, we will deal with spinless fermions described 
by annihilation and creation operators, Ci and ct, re¬ 
spectively. In this case too, we are interested in two- 
point functions that characterize the correlations and the 
causal response of the model. With two (annihilation and 
creation) operators at our disposal, the correlation func¬ 
tion becomes a 2 x 2 matrix defined as 

j)]^p = (Tc,ait)c]p{0)'j , (3) 

where a, (3 G { 1 , 2 } with (cn ca) = {ci c]) and the 
fermionic operators given in the Heisenberg picture—the 
long-ranged (1/|* — lattice Hamiltonian is speci¬ 

fied in Sec. IV. T is the time-ordering operator defined as 
rO'{t)O{0) = O'{t)O{0) and rO( 0 )O'(t) = -O'{t)O{0) 
for t > 0 and fermionic operators O and O' (site indices 
suppressed). Similarly, the response function is defined 
as 

i [D„(<, i - j)]^p = 0(t) ^ [c^a{t), c]^(0)]+^ , (4) 

where the brackets with the subscript -I- denote the anti¬ 
commutator. Here too, a Lieb-Robinson bound dictates 
that the response function is exponentially suppressed 
outside a linear light cone for short-range interactions 
[ 21 ]. 


III. Transverse-Field Ising Model: Scalar Field 

In this section, we consider the critical and causal prop¬ 
erties of the transverse-field Ising model with long-range 
interactions. Let us first consider the TFIM with nearest- 
neighbor interactions in one dimension, 

H = -Y.{StS^+,+gS^), (5) 

i 

with the Pauli operators. This model undergoes 

a quantum phase transition at g = 1 from an ordered 
phase {g < 1 ), where the Z 2 symmetry of the model is 
broken and {Sf) 7 ^ 0, to a disordered phase {g > 1) where 
(S'?) = 0. Near the critical point, the long-wavelength 
behavior of the model can be described by a field theory 
in the continuum [34], with the Euclidean action 

I = J dr J dx {dr(j)Y + {9x4>Y + ■ (6) 

Here we have rescaled time and spatial coordinates to 
normalize the coefficients of the derivative terms, defined 
g g — 1 characterizing the vicinity to the critical point. 


and u as the interaction strength. The above action de¬ 
fines the (j)'^ field theory in two-dimensional Euclidean 
space, with (j) a coarse-grained field representing (S^); 
both are the Ising order parameters which measure the 
broken symmetry across the phase transition. Indeed, 
one can see that the scaling dimension of S*, i.e. how 
its correlations behave under rescaling of space and time 
coordinates, is consistent with that of the Wilson-Fisher 
fixed point in two dimensions and with N = 1 component 
[34]. 

Next, we consider the transverse-field Ising model with 
long-range power-law interactions. 




SfS^ 








( 7 ) 


where the exponent is defined in terms of cr > 0 ; the 
latter is assumed to ensure a well-behaved thermody¬ 
namic limit, in which the energy remains extensive as the 
system size approaches infinity. (The anti-ferromagnetic 
version of this model, studied in Ref. [35], is not identi¬ 
cal to the ferromagnetic model owing to the long-range 
interactions). A numerical study of the phase diagram of 
the above model is performed in Ref. [36]. We also note 
that long-range interacting classical Ising models (with¬ 
out the transverse magnetic field but at finite tempera¬ 
ture) have also been studied extensively in one dimen¬ 
sion [37, 38]. The model in Eq. (7) undergoes a quantum 
phase transition similar to the short-range TFIM, but at 
a cr-dependent critical coupling strength g, whose pre¬ 
cise value will not be important for our purposes; we are 
rather interested in the universal properties of this model, 
namely the scaling of various correlation functions with 
distance and time, while the precise value of the critical 
g can only affect constants of proportionality. The cor¬ 
responding universality class also includes models with 
additional ferromagnetic nearest-neighbor, next-nearest- 
neighbor, or higher but finite range terms, added to the 
Hamiltonian in Eq. (7). 

Mapping onto the continuum generates long-range in¬ 
teractions in (j) and its powers provided the original sym¬ 
metries of the Hamiltonian are respected. For example, 
(j){T,x)(l){T,y)/\x - and (l){T,x)(j)^{T,y)/\x - y\^+'^ 

are odd with respect to 4>{t,x) —>■ —4>{t,x), i.e. a spin 
flip —Sf, similar to the first term in but 

(/)(t, x)(j)'^{T, y)/\x — y\^~^'^ is not. For cr > 0, all long-range 
terms beyond the quadratic order in (j) and all terms in¬ 
volving spatial derivatives beyond the (local) quadratic 
gradient term are irrelevant in the sense of RG; additional 
insertions of the field will make the corresponding term 
less relevant in the sense of RG, hence (j){x)(l){y)^/\x — 
is less relevant than (f)[x)(j){y) /\x — , for ex¬ 

ample. The full action then includes a quadratic piece, 
containing the long-range term ^(r, x)(j>{T, y)/\x — y\^'^'^^ 
plus the local (j)^ interaction. The former can be more 
conveniently cast in Fourier space (in imaginary time) 
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[39], such that 


A. Quadratic model 


I = + J dr J dx u(j)'^, ( 8 ) 

=Jduj jdq {uj'^ + g + g)p . (9) 

The parameter g now defines the distance from the 
critical point of Eq. (7); here, we will restrict our¬ 
selves to p > 0 , that is, we consider either the criti¬ 
cal point, or the paramagnetic side of the Ising transi¬ 
tion. The dispersion relation at the quadratic order is 
simply uJa{q) = \/g + q'^ + i?cr|q|'^ ■ Note that the dis¬ 
persion relation computed by Fourier transforming the 
coupling constants of the Ising term in Hg. is more com¬ 
plicated (leading to an expression in terms of polylog¬ 
arithms). Here, in the interest of describing the long- 
wavelength physics that the continuum description is 
suited to, we have only kept the leading low-g analyti¬ 
cal (q^) and non-analytical (jqj'^) terms (non-analytical 
terms with a power smaller than a cannot appear since, 
upon the inverse Fourier transform, they would give rise 
to power-laws that decay slower than the original power- 
law interaction ^ l/r^+'^); both will play a crucial role in 
the following sections. Higher orders of momenta would 
create finer, or faster-decaying, features, and are thus ig¬ 
nored. 

The two-point functions in Eqs. (1) and (2) can be 
described in the continuum as 

iGa{t,x) oc {T(j){t,x)(j){0,0)), (10) 

iD„{t,x) oc 0(<)([(^(t,x),(^(O,O)]), (11) 

where (j) denotes the operator-valued field in the Heisen¬ 
berg picture. In these equations, we have only 
given the proportionality relations, as we are ulti¬ 
mately interested in scaling relations but not the precise 
coefficients—the coarse-graining of the lattice operators 
to their continuum counterparts will introduce nonuni- 
versal lattice-constant-dependent coefficients, while the 
long-distance/time behavior are unaffected apart from 
overall coefficients. We note that, although Lieb- 
Robinson bounds do not apply directly to the contin¬ 
uum model being studied, sufficiently close to criticality 
the continuum response function gives a quantitative de¬ 
scription of the lattice response function on length scales 
much larger than the lattice spacing, and hence encodes 
the correct asymptotic scaling with space and time of the 
response function for the underlying lattice model. 

In the remainder of this section, we study the con¬ 
tinuum description of the Ising model with long-range 
interactions at or away from criticality within the RG 
approach. In all cases, we restrict ourselves to the vicin¬ 
ity of the critical point, where the correlation length is 
large compared to the lattice spacing, and thus a contin¬ 
uum description and the RG treatment are well justified. 


To gain some insight into the causal structure and cor¬ 
relations of the long-range interacting model, we first 
study the quadratic part of the action, but consider the 
effects of the interaction via RG in the next section. At 
the quadratic order, the two-point correlators above com¬ 
pletely characterize the system, and can be directly ex¬ 
tracted from the dispersion relation as (see App. A) 

Ga{R) oc f cos{qR) (12) 

Jo 

for the correlations at equal times [G„{R) = Gcr( 0 ,i?)] 
and 

pOO 1 

Da{t,R)(x / —^ sm{oJa{q)t)cos{qR) (13) 

Jo ^cr[q) 

for the response function. Here we denote the distance 
by i? = I a; I, and have only given the proportionality rela¬ 
tions, as we are ultimately interested in scaling relations 
and not the precise coefficients. 

We start by making a simple observation that, for 
a < 2, the long-range term (jgj'^) appears to be more 
relevant than q^, in which case one can simply drop the 
latter from Eq. (9). However, a more careful treatment 
reveals that, at least at criticality, the long-range inter¬ 
action, considered as a perturbation on top of the short- 
range interacting model in Eq. (6), can be dropped for 
values of (T > 2 — t^sr = 7/4 where ? 7 sr = 1/4 is the 
anomalous dimension of the field (j) in the short-range 
interacting model ( 6 ). In the remainder of this subsec¬ 
tion (at the level of the quadratic action), we shall ig¬ 
nore this complication, but discuss it in some detail in 
Sec. HIB. We also remark that a specific feature of the 
model in Eq. (7) is that 5^ > 0 for 0 < tr < 2, and 
Ho- < 0 for 2 < cr < 4, a pattern repeated periodically; 
this can be easily seen by computing the Fourier trans¬ 
form of the \/\i — power law in Eq. (7). The fact 
that Ba is positive when it is relevant is indeed assur¬ 
ing as the Hamiltonian is bounded from below. On the 
other hand, when the corresponding term is not relevant, 
we may have B^j < 0 which does not pose a problem as 
one must impose a high-momentum cutoff naturally pro¬ 
vided in lattice systems. Below, we study these two cases 
separately. 

For the sake of comparison, we first quote the two- 
point functions in the absence of long-range interactions 
at or away from criticality. First, the correlation func¬ 
tion for the quadratic model at criticality is given by 
Gsk{R) logR, while the full (short-range) interacting 
model yields Gsr(R) [34]. Away from critical¬ 

ity, the correlation function decays exponentially beyond 
a length scale defined as the correlation length i.e., 
Gsr(7?) ~ exp(—i?/^). The response function is given 





6 


by 


DsR{t,R) 


Q{t — R), 0 = 0, 

0(t)lm [i^o i^/eVR^ - i2)] , g>0, 


( 14 ) 

that is, it vanishes identically outside the linear light cone 
t = R, and, for ^3 = 0, it is simply a constant within 
the light cone. We note that the response functions in 
Eq. (14) would be modified if interactions were included, 
and thus serve only as approximations to the response 
functions for the short-range interacting model. 


1. 0 < a < 2 at criticality 

At the critical point, the dispersion relation can be ap¬ 
proximated as uia{q) « where we have set B^r = 1. 

The correlation function can be computed from Eq. (12) 
via a simple rescaling of g by 1/i?, yielding 

( 15 ) 

The exponent of this power-law defines the scaling di¬ 
mension of the field at the quadratic order; however, it 
will receive corrections in the course of RG, as we shall 
discuss in Sec. IIIB. 

The response function is obtained from Eq. (13) via a 
similar rescaling, yielding 



R 


Figure 1: Contour plot of the (unnormalized) response func¬ 
tion D 2 / 3 {t, R). The response function is computed via the 
scaling function 02/3 • The causal region exhibits a nonlin¬ 
ear behavior with a finger-like pattern extending beyond the 
linear light-cone structure of short-range interactions. The 
numerical evaluation of the response function with the lat¬ 
tice dispersion relation and momentum cutoff produces an al¬ 
most identical plot. The (red) solid line represents the curve 
t ~ with cr = 2/3. 


first local maximum in time around 


Dcr{t,R) 


1 r dq 

Jo 

1 


qal2 

t 

R '^/2 


sin 



t 

R '^/2 


cosq 

(16) 


where we have defined the scaling function in the last 
equality. Note that, by extending the momentum integra¬ 
tion to infinity, we have derived a simple scaling relation 
for D„{t, R). We point out that a similar scaling relation 
also emerges for the time-ordered correlation function 
Ga{t,R), though with a different functional form. Nev¬ 
ertheless, one should keep in mind that the above expres¬ 
sion applies only to regions well outside the linear light 
cone R = t since we have dropped the momentum cutoff 
as well as the high-momentum modes (including q^) from 
the dispersion relation. We also note that the multiplica¬ 
tive power-law in Eq. (16) is identical to Eq. (15) simply 
due to the scaling dimension of the field (/>; apart from this 
power-law required for dimensional reasons, the informa¬ 
tion about the causal behavior is encoded in the scaling 
function g^- For small arguments, ga{s) indepen¬ 
dent of (T, which results in Dcr{t,R) ~ Thus, 

at fixed t, the response function decays as 1/1?^+°' with 
distance, consistent with the demands of the Hastings- 
Koma bound [22] applied to Her [74]. The function ga{s) 
increases monotonically with s up to s ^ 1, but, beyond 
this point, exhibits an oscillatory behavior to be further 
discussed below. Hence, the response function at a dis¬ 
tance R away from a local quench at t = 0 reaches the 


t - . (17) 

This implies that the causal region obeys a power law, 
in sharp contrast with short-range interactions. A nu¬ 
merical evaluation (at the quadratic order) with the 
full dispersion relation shows that the response function 
branches outside the linear light cone polynomially as 
t ~ see Fig. 1. In fact, by setting Dcr{t,R) to a 

small constant, we find polynomial fingers reaching out¬ 
side the linear light cone t = R, which are, however, finite 
in extent as the signal weakens while it propagates. In 
general, the effect of a local quench decays at long times, 
and thus contours of the response function—due to a lo¬ 
cal perturbation—are typically finite in extent. The fact 
that the response function for short-range interactions 
[Eq. (14)] does not decay at long times is an artificial 
feature of the quadratic action and would not be 
the case in the presence of interactions, in higher dimen¬ 
sions, or sufficiently far from the critical point (where RG 
irrelevant terms that break relativistic invariance must be 
included). Therefore, to correctly identify the causal be¬ 
havior, we disregard the decaying power-law function out 
in front (in this case, 1/7?^“°’^^), and focus on the scaling 
function g^. 

The power law in Eq. (17) identifies the dynamic crit¬ 
ical exponent that characterizes the relative scaling of 
time with respect to space. For cr > 2/3, RG calcu¬ 
lations produce corrections to the mean field value (i.e. 
calculated at the quadratic order) of the exponent above. 
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as we shall discuss in Sec. IIIB. Scaling relations of the 
form (16) are generic beyond the quadratic model. How¬ 
ever, within the quadratic model, an explicit form for the 
scaling function gcr can be obtained: For all 0 < ct < 2 it 
undergoes periodic oscillations in units of of 
with an asymptotic envelope function 

Amp[g<,(s)] - s(i—)/(2—) . ( 18 ) 

For example, for ct = 1, we have g{s) ~ cos(s^/4) -|- 
0{s~^). Compared with numerics, the scaling function 
g^r is quantitatively accurate for the low lying fingers, but 
also captures the qualitative features of higher ones. 

2. a > 2 at criticality 

In this regime, long-range interactions are not relevant, 
and one thus expects to recover the same asymptotic form 
of the correlation function as well as the linear light-cone 
structure of short-range interactions. However, the slight 
deviation from a linear light cone can be quantified, and 
shown to take a universal form. For convenience, we re¬ 
strict ourselves to 2 < ct < 4, since in this range higher- 
order analytical terms in momentum (g^, g®, • • •) can be 
dropped compared to |g|'^. The physics at long times and 
distances is dominated by its behavior at low momen¬ 
tum, where the dispersion relation can be approximated 
as a;cr(g) « |g| — Igl*^”^. Here we have dropped the coeffi¬ 
cient of the non-analytic term, but kept the correct (neg¬ 
ative) sign inherited from B^. We also approximate the 
frequency in the denominator of Eq. (12) by uja-{q) « g. 
This set of approximations is valid in the vicinity of the 
linear light cone t = R, which is the main focus here. For 
space-time points well inside or outside the light cone, 
the response function is generally not captured by the 
following scaling relations. With the above approxima¬ 
tions, the response function can be recast as 



Figure 2: (a) The scaling function /o- characterizing the re¬ 
sponse function as a function of s = (t — for 

a = 5/2. The function /o- approaches a constant inside the 
linear light cone (s > 0), but decays rapidly outside the light 
cone (s < 0). (b) The first local maximum of the response 
function (solid curve) is bent inside the light cone (dashed 
line). / characterizes the deviation from the linear light cone. 
For long times, ^/t —> 0, and the causal region becomes a 
sharp linear cone. 

This implies that the response function approaches a con¬ 
stant (zero) inside (outside) the linear light cone. The 
crossover between the two limits only depends on R and 
t in the particular combination in Eq. (19). To find the 
causal region, we set the response function to a (small) 
constant which then yields a constant value for the argu¬ 
ment of the function f^, 

t — R 

Eor sufficiently long times, we have <^t — R near 

R Ki t, which yields a nearly-linear light cone; however, 
the causal region is not a sharp cone as it spreads out over 
a region of size <c t. The deviation of the causal 

region from the linear light cone can be characterized by 
^ = t — R-, the above considerations yield 


Da{t,R) 


‘ [ — { sin[g(t - R)-q'^ 4] 

Jo 9 ^ 

-I- sin[g(t -I- i?) - 


=h 


t — R 

^l/(cr-l) 


+ I<y 


t + R 


At sufficiently long times {t ^ 1), the argument of the 
last term is large since <C t for ct > 2. In this 

limit, one can see that /o-(s) —>■ const as s —)■ -l-oo. There¬ 
fore, the response function can be written as 


At long times, ^/t ^ 0 for ct > 2 , and the linear light 
cone becomes exact. For the quadratic model considered 
in this section, one can sometimes (for certain values of ct) 
obtain an analytical expression for the scaling function in 
Eq. (19). We find that the causal region is slightly bent 
inside the linear light cone [with ^{t) ^ gee 

Eig. 2. The fall-off outside the light cone is monotonic, 
with asymptotics governed by 


Da{t,R) 


fa 


t — R 


(19) 


fa{s)^\s\^-^. ( 22 ) 


where the function fa-{s) = fa{s) + fai+oo) has the prop¬ 
erties 


fa{s) 


const, s —>■ -boo , 
0 , s —>■ —oo . 


( 20 ) 


3. All a > 0 away from criticality 

We first compute the correlation function Ga{R) away 
from criticality for all ct > 0 not being an even inte¬ 
ger. With short-range interactions, this function decays 




















exponentially beyond the correlation length; however, 
long-range interactions always induce a power-law de¬ 
cay of the correlation function. To be completely gen¬ 
eral , we can consider a dispersion relation of the form 
Wcr(g) = \/l + F{q^) -I- Ba\q\'^ where F{q^) is a polyno¬ 
mial function of [with F{q^ —>■ 0) = 0] such that the 
expression under the square root is nonnegative; while 
we usually truncate the function F{q^) at the quadratic 
order to capture the long-wavelength physics, it gener¬ 
ally has a more complicated, but analytic, form. Equa¬ 
tion (12) can be cast as an integral over q S (— 00 , 00 ) 
with the substitution cos(gi?) One can then de¬ 

form the contour of integration to the upper-half plane in 
the complex plane. However, the function \qY’ should be 
treated separately as an analytic function for q ^ 0. We 
thus have to find the zeros, or branch points, of uJa-{q) in 
the upper half-plane. For a value of cr that is not an even 
integer, there are no zeros on the real or imaginary axis. 
The branch points are denoted by q*, which, with a con¬ 
venient normalization, can be assumed to be of the order 
of unity, and specifically Img* ~ 1. The contribution of 
a branch cut is exponentially suppressed as 

X(i?)e-(l”«*)«, (23) 


up to a multiplicative polynomially decreasing function 
x(i?). For large R, the branch cut only contributes near 
q ^ q*, thus ujaiq) ~ V? “ 9* and x(i?) ~ / dq^e~^^ ^ 

Ijy/R. On the other hand, the integral along the imagi¬ 
nary axis gives 


G.(i?) 



dq Im , 

y/l + F{-q^) + 


dq 




1 


(24) 


up to exponentially decaying corrections in Eq. (23). De¬ 
pending on the relative ratio of the coefficients of the ex¬ 
ponential and power-law terms in Eqs. (23) and (24), we 
may find an exponentially decaying correlation at inter¬ 
mediate length scales followed by a power law at large 
distances, consistent with Refs. [25, 40, 41]. In the last 
step, we have expanded the denominator around g = 0 to 
obtain the large-distance behavior, which yields the exact 
asymptotic form of the correlation function. Replacing 
q'^ by an analytic expression such as g” with n an even 
integer, this argument breaks down as lTa(l/L>jq=i\q\) = 0 
on the imaginary axis. 

The exponent characterizing the decay of the correla¬ 
tion function away from criticality [Eq. (24)] thus coin¬ 
cides with the 1 -|- cr exponent of the long-range interac¬ 
tions. While this result is derived at the quadratic level 
of the action, it holds more generally even in the pres¬ 
ence of interaction terms, as we shall discuss in the next 
section. 

Next we study the response function. Away from crit¬ 
icality, scaling relations take a more complicated form, 
and analytical expressions are less available. Instead we 



Figure 3: The norm of the (unnormalized) response function 
\Di{t,R)\ away from criticality evaluated with the full dis¬ 
persion relation uJi{q) = for q £ (0, tt) (with 

= 1 and the g-dependence obtained from the polylog func¬ 
tion of the lattice dispersion relation at a = 1). The analytical 
expression in the text produces a very similar plot sufficiently 
away from the vertical line {R — 0). The response function 
exhibits a linear causal behavior. 


shall resort to a simple analysis of the group velocity 
Vq = duJcr{q)/dq. A semiclassical approximation sug¬ 
gests that the causal region should be linear when the 
velocity is bounded from above [42], see also [26]. The 
group velocity at high momenta is typically bounded due 
to the short-wavelength cutoff; however, for sufficiently 
long-range interactions, it may diverge as g —>■ 0. At 
criticality, the semiclassical picture correctly predicts a 
linear causal behavior for cr > 2. Away from critical¬ 
ity, the dispersion relation is ^^(g) = a/I + g^ -I- B„\q\'^, 
where we have set g = 1. An inspection of the group 
velocity then shows that, in the limit where g —f 0, it di¬ 
verges for CT < 1, but approaches a constant (or zero) for 
CT > 1. For the special value of cr = 1 at the borderline 
between the linear and non-linear causal behaviors, the 
dispersion relation is approximately Wi(g) « F 2|g|, 
where we have set i?o- = 2 for notational convenience. In 
this case, one can find an exact analytical expression for 
the response function as 


Di{t, R) ^ Re-^ 


6 -%(R^+R)/2R 

Vr 


X 



R-t\ 


{t -f -t) 


(25) 


where erf is the error function. The expression inside the 
curly brackets falls sharply beyond a straight line t = R, 
which thus indicates a linear boundary of the causal re¬ 
gion. The response function, i.e., the real part of the 
expression in curly brackets, is highly oscillatory, but ex¬ 
hibits the same causal behavior as the expression inside 
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|Da (R“,R)| 
1 



- a 

- a 

. a 


Figure 4: The norm of the (unnormalized) response function 
\Da{t, -R)| evaluated at t for cr = 1/3 away from critical¬ 
ity. The three plots correspond to different choices of a = cr/2, 
a, and 1 (with a = 1/3). The dispersion relation is taken 

as uji/siq) = (^) (^) 1 ^ (O.^^) 

(with g = 1 and the g-dependence obtained from the polylog 
function of the lattice dispersion relation at a = 1/3). For 
the exponent a = ( 7 / 2 , we find a simple power-law decay of 
the response function with R, suggesting 2 : = a 12. 


the brackets. Furthermore, the amplitude of the response 
function decays as 1 /E? for fixed t outside the linear light 
cone, consistent with the Hastings-Koma bound [22] for 
the lattice model (7). The response function computed 
with the full dispersion relation, plotted in Fig. 3, clearly 
exhibits a causal region with a linear boundary. We also 
note that response functions are generically highly os¬ 
cillatory away from criticality even for the short-range 
interacting model [Eq. (14)]. 

Next we consider cr < 1. In the absence of analytical 
expressions, it is more difficult to fully characterize the 
causal behavior. Naively looking at the response function 
could be misleading as setting Dc,{t, R) =const gives rise 
to a vertical line at long times in some cases, and does 
not clearly identify the causal behavior. Instead we set 
t ^ i?“ in D„ {t, R) and plot it as a function of R. We can 
then identify 2 as the value of a for which the response 
function obtains a simple power-law form. The function 
Dcr{R°‘, R) is plotted in Fig. 4 for cr = 1/3 and three dif¬ 
ferent values of a. Clearly, only forQ; = cr/2 = l /6 does 
this function have a particularly simple form, falling off 
with the distance as a power-law l/i?^+°’ for ct = 1/3. As 
expected, the associated exponent is identical to the one 
characterizing the decay of Ga{R) at long distances in 
Eq. (24). These observations are consistent with the scal¬ 
ing form Dcr{t,R) ^ with a = 1/3 

and 2a- a scaling function, at least for the set of pa¬ 
rameters studied here. The scaling function 2a defines 
the relative scaling of space and time coordinates, giv¬ 
ing a dynamic critical exponent z — for values of 
0 < cr < 1; and in this sense the causal behavior is non¬ 
linear. This is identical to the dynamic critical exponent 
of Sec. IIIAl computed for 0 < cr < 2 at criticality, 
while, in the noncritical case studied in this section, the 
above identification is valid only for 0 < cr < 1. We stress 


that this conclusion is based on our numerical results for 
a limited range of parameters; an extensive investigation 
of this regime should be carried out to confirm its uni¬ 
versal scaling properties. 

--all 

--0 

^ B. Beyond the quadratic model: A 

renormalization-group study 

To study the full interacting model in Eq. ( 8 ) , we 
first recall that long-range interactions are irrelevant for 
a > 2 — ? 7 sr = 7/4, in which case the renormalization 
group and the universal properties and exponents at crit¬ 
icality are those of the usual model in one higher di¬ 
mension than the original model. The borderline value 
of cr = 2 — rysR arises due to the following consideration. 
We can treat the long-range interaction jqj'^j0(a;, q)\‘^ as a 
perturbation on top of the short-range interacting model 
in Eq. ( 6 ), for which the scaling dimension of the field is 
[■ 5 ^] = »?sr/ 2. One can then easily see that the scaling di¬ 
mension of the long-range interacting term becomes neg¬ 
ative for values of ct > 2 —rysR- This feature also arises in 
classical long-range interacting systems [43] , see also [44] , 
and also [45] for a recent review. We thus first write the 
action for 0 < cr < 2 — rygu in real space and (imaginary) 
time in d spatial dimensions as 

I = [a {drcjdf' + 

+ ff ^ ’ 

with the (possibly normalized) coefficients A, Bn, and g. 
As a first step, we rescale the coordinates and the field 
as 

X —)• x' = x/5, T ^ t' = , (j) ^ (j)' = b°'4' j 

where b is an arbitrary constant greater than unity; note 
that the time and spatial coordinates are scaled differ¬ 
ently. The dynamic critical exponent, z, and the scaling 
dimension of the field, a, are to be determined. Sim¬ 
ple dimensional analysis yields the scaling dimension of 
various terms in the action as 

[A] = —z + d — 2a, [Bn] = z — a + d — 2a, 

[p] = z + d — 2a, [it] = z + d — Aa . (26) 

To find the exponents at the level of mean field, we set 
[A] = [Bn] = 0, that is, we require the quadratic part of 
the action (with the exception of the ‘mass term’ g) to 
be invariant under rescaling. We then obtain 


CT 

"= 2 ’ 


^ = 1 - 2 - 


(27) 


where we have defined the anomalous dimension 77 of the 
field via a = (d — 1 + r])/2. Indeed, with d = 1, these 
are the same exponents that describe the two-point func¬ 
tions in Sec. IIIAl. To go beyond mean field, we shall 
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resort to the so-called epsilon expansion. To this end, we 
first derive the upper critical dimension beyond which the 
interaction becomes irrelevant, i.e. [u] < 0. Demand¬ 
ing that [u] — Q at d = du, one finds du = 3tT/2 [.39]. 
Therefore, for our one-dimensional model, the interac¬ 
tion term can be dropped when a < |, in which case the 
field-theoretical arguments indicate that the full model 
is faithfully represented by the quadratic part, and the 
mean-held exponents become exact. On the other hand, 
for cr > 2/3, these exponents may be modihed due to 
fluctuations, and corrections to them can be obtained as 
an expansion in e = — d = 3 (t/2 — 1 [46, 47]. 

Before proceeding, we remark that the RG procedure 
produces only analytical terms, and thus cannot renor¬ 
malize non-analytical terms in the action. Therefore, the 
coefficient Bu in the action is not renormalized, which 
leads to an exact non-renormalization condition from 
Eq. (26) by setting [Ro-j = 0 (cast in terms of the anoma¬ 
lous exponent rf) 

r] = 1 + z — a , (28) 

valid to all orders of perturbation theory. Thus Hnding 
one of the two exponents (77 or z) completely determines 
the other one. (For long-range interacting classical sys¬ 
tems, i.e., in the absence of the imaginary-time direction, 
the exponent 77 assumes the hxed value 77 = 2 — ct [48].) 

The e-expansion can be organized as a series of loop 
diagrams. The two-point function is renormalized by the 
loop diagrams depicted in Fig. 5 to the two-loop order. 
However, the one-loop diagram does not contribute to the 



Figure 5: The one- and two-loop diagrams contributing to the 
two-point function. 

exponents 77 or z since the loop correction is independent 
of the (imaginary) frequency and momentum of the ex¬ 
ternal lines. Therefore, we should consider the two-loop 
diagram, which is significantly harder to compute. We 
leave the details of RG to App. B, but quote the result of 
the second-order e-expansion for the exponents z and 77 . 
Defining Az and A 77 as the difference of the exponents z 
and 77 from their mean-field values in Eq. (27), we find 

A 77 = Az = (r(o')e^-I-C>(e^), (29) 

where <;(a) is a rather complicated expression reported 
in App. B, but is approximately reproduced by c(o’) ~ 
1/[24(1 -I-cr^)]. (For the special case of tr = 1, this is con¬ 
sistent with the result obtained in Ref. [49], see App. B.) 
Gomputing the critical exponents in d = 1 dimension, we 
have e = 3 (t/ 2 — 1 as explained above. We briefly note 
that, while e can be of order 1 , e—expansion has been 
remarkably successful even for e = 1,2 [50]. The expo- 


n z 



Figure 6 : The critical exponents 77 and z as a function of a. 
These exponents assume their short-range values 77 = 77 SR = 
1/4 and z = 1 for a > 7/4. For a < 2/3, they stick to 
their mean-field values 77 = 1 — a/2 and z = cr/2, and, for 
2/3 < a < 7/4, their values are computed up to the two- 
loop order (the dashed lines show the mean-field exponents). 
The inset shows Az = z — a/2 and A 77 = 77 — (1 — a/2), 
the difference of the critical exponents from their mean-field 
values, up to the two-loop order. 


nents 77 and z as a function of a are plotted in Fig. 6 . 
For a < 2/3, field theory indicates that these exponents 
stick to their mean-field values, while for a > 7/4, they 
assume their short-range values. For intermediate val¬ 
ues of 2/3 < a < 7/4, the two-loop correction gives a 
slight deviation from the mean-field prediction (the lat¬ 
ter shown as the dashed line in the range 2/3 < tr < 7/4 
in Fig. 6 ). The two-loop corrections to the exponents 77 
and z are plotted in the inset. Importantly, there appears 
to be a discontinuity in the values of both exponents at 
a = 7/4; however, it may very well be the case that 
higher-order corrections in e-expansion make the jump 
disappear. While the two-loop correction is rather small 
at the transition near cr = 7/4, it gives corrections to both 
exponents towards their short-range values, and specif¬ 
ically brings the model closer to the ‘relativistic’ point 
where z = 1, and, by virtue of Eq. (28), also 77 = 1/4. 
In the context of long-range interacting classical systems 
(in the absence of the imaginary-time direction), it was 
shown that there is no such discontinuity [43], but this 
result has been the subject of further scrutiny recently, 
see Refs. [26, 51, 52]. Further analytical and numerical 
work in the quantum context should also be worthwhile. 

Away from criticality, we note that RG does not mod¬ 
ify power laws obtained at the quadratic order. This 
indicates that the exponent in Ga{R) ^ 1 /i?^^'’’, com¬ 
puted for the noncritical model at the quadratic order, 
does not receive any corrections from RG. 
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IV. Fermions 


In this section, we study the critical and causal proper¬ 
ties of a model of spinless fermions on a chain with long- 
range hopping and pairing terms. To set up the problem, 
we first introduce a short-range quadratic lattice model 


H = - '^{c\c^+l + c\^^Ci) + J{c\c\^^+c^+lC^)-^lc\ci, 

(30) 

with nearest-neighbor hopping and pairing (with the am¬ 
plitudes 1 and J, respectively), and a chemical potential 
/i. The long-wavelength physics of this model is captured 
by a continuum field theory with the Euclidean action 
[34] 


(31) 

where 4' is a Grassman field. The spatial coordinate 
is rescaled to normalize the coefficient of the gradient 
term, and the constant A is generically related to the 
parameters in the lattice model; for example, for J = 1, 
we have A = 2 — /r. We shall take the continuum de¬ 
scription as our starting point, and are not concerned 
with the relationship of its parameters to those in the 
lattice model. In principle, the Hamiltonian in Eq. (30) 


may also include interaction terms such as c|cic|^_|_]^Ci+i. 
Mapping to the continuum, such an operator induces an 
interaction term that must involve spatial gradients (as 
a result of the Pauli exclusion principle, = 0). From 
the scaling dimension of the fermionic held, ['['] = 1/2, 
one can easily check that all allowable interaction terms 
are irrelevant in the sense of RG [34], and thus Eq. (31) 
is the exact (in the RG sense) continuum description of 
the lattice model even with interactions. We also remark 
that the fermionic held theory in Eq. (31) is relativistic 
in the sense that the dynamic exponent is z = 1, indi¬ 
cating a linear causal region. The action can be cast in 
Fourier space as (integral over momentum and frequency 
suppressed) 




—iu! + A iq 
—iq —ioj — A 



where w is the Fourier variable corresponding to the 
imaginary time. Setting the determinant of the above 
matrix to zero, we hnd the dispersion relation in imagi¬ 
nary frequency w = iuj{q) with uj{q) = ■>/A^ -|- q'^. 

We are interested in the consequences of adding long- 
range terms to a fermionic model in Eq. (30). For exam¬ 
ple, one can add both long-range hopping and long-range 
pairing terms. 

Hopping : 

Pairing : 


1 

l*-/l 
1 


X/ I,- _ tII+o- ’ 
i¥=3 

E 


Q 


i>j 




Again, we shall assume a > 0—for — 1 < cr < 0, our treat¬ 
ment requires special care, and will be extended in future 
work. Long-range interactions, beyond the quadratic or¬ 
der in fermionic operators, can also be considered, though 
they do not affect the possible RG-relevant contributions 
to the continuum theory. Mapping these expressions to 
the continuum, the matrix in Eq. (32) becomes 

f-iuj + A- q'^+Ha(q) iq + Vaiiq) \ /gg-, 

y -iq + Vai-iq) -iuj - A + q'^ - Ha{q) j ' 

in which 'Hcr(g) ^ \q\'^ and Va{iq) ~ (with 

± corresponding to q ^ 0 ) denote low-momentum ex¬ 
pansion of the long-range hopping and pairing terms, re¬ 
spectively. In the case of long-range hopping, 'Pa{iq) 

which also diverges for odd-integer val¬ 
ues of cr. In this case, we shall restrict ourselves to 
non-integer tr > 0. In the above equation, we have 
also included a quadratic term in q (with a normal¬ 
ized coefficient) on the diagonal, even though it is ir¬ 
relevant compared to the linear-momentum term; the 
former will be important in some cases at criticality 
(A = 0). The dispersion relation is now given by 
Wa(g) = ^/[A - q^ +na{qW‘ + \n + 'Pc{iq)\^. 

Comparing the exponents of the momentum- 
dependent terms in Eq. (33), we find that long-range 
terms (~ l?]*^) are either relevant or irrelevant compared 
to the linear-momentum term (~ q) depending on 
whether cr < 1 or cr > 1, respectively. Therefore, the 
dynamic critical exponent is z = cr when 0 < cr < 1 , 
and z = 1 for cr > 1. Furthermore, the scaling dimen¬ 
sion of the fermionic field is fixed by the linear time 
derivative [the hrst term in Eq. (31)] and remains the 
same, [dr] = 1/2, independent of cr. Hence, by the 
same arguments that we applied to the local model, all 
possible long-range non-linear terms are irrelevant. The 
scaling dimension of the Fermi field is discussed in more 
detail in Sec. IV 1. 

Next we define the two-point functions of the Fermi 
field in the continuum in simple analogy to their lattice 
definitions, Eqs. (3) and (4). The correlation function is 
given by 

i[Gcr{t,x)]^p cx {|T^a(t,a;)^^( 0 , 0 )^ , 

where ^ and are the operator-valued Fermi helds 
in the Heisenberg picture. We recall the definitions 
'I'^^ and cr,/3 G {1,2}. Similarly, the 
response function is defined as 

i [D„(t, x)]^p oc 0 (t) [4'a(t, x), 4' 1 ,( 0 , 0 )]+^ . 

In the above equations, we have only given the propor¬ 
tionality relations, as we are only interested in scaling 
relations and not the (nonuniversal) coefficients. 

With the fermionic action being quadratic, the two- 
point functions can be extracted from the dispersion re¬ 
lation. The correlation function at equal times, G„{R) = 
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Ga{t —)• 0, R) (with R = |a:| the distance between the two 
points), can be written as [75] 

G.(i?) cx ^ 2 ^ G,iR) + GAR) ■ (34) 

Here GaiR) is the correlation function defined in 
Eq. (12), but with uj^iq) now being the dispersion of 
the fermionic model, and the dispersion ujuiq) of the 
fermionic model, and 


GAR) oc 



RAq) 

Rai-iq) 


R-yiiq) 

-RAA 


(35) 

As before, we have given the proportionality relations as 
we are interested only in scaling relations and not the 
precise coefficients [the factor of 1/2 in Eq. (35) is rela¬ 
tive to the first term on the right-hand side of Eq. (34)]. 
In the above equations, the diagonal matrix elements give 
the normal correlation function (('I'4'l) and its complex 
conjugate), while the off-diagonal matrix elements define 
the anomalous correlation function ((4''P) and its com¬ 
plex conjugate). Similarly, the response function can be 
written as 


Do.(t, R) oc 



Or 

idt-A- d\ 


Da{t, R)+T)A^t, R) 


(36) 

where Dcr{t,R) is the response function defined in 
Eq. (13), but with uj„{q) now being the dispersion of the 
fermionic model, and to the dispersion relation oJa^q) of 
the fermionic model, and 


indicating that the response to an infinitesimal pertur¬ 
bation travels at the speed of ‘light’. Note that we have 
dropped the second spatial derivatives; this and higher- 
order derivatives may be included, but will only slightly 
smear out the delta function. Away from the critical 
point, the response function remains zero outside the 
light cone, but inside the light cone it has a compli¬ 
cated form given by the first term in Eq. (36), where 
DsR{t,R) should be substituted from Eq. (14) at > 0 
with Aff A. In what follows, we undertake a de¬ 
tailed study of the fermionic field theory with long-range 
hopping and pairing, which is rather rich despite the ab¬ 
sence of interactions. We will explore the linear/non- 
linear causal behavior, expressing it in terms of explicit 
scaling relations whenever possible. In all the cases stud¬ 
ied below, we shall restrict ourselves to the vicinity of 
the critical point, where the correlation length is large 
compared to the lattice spacing, hence the validity of a 
continuum description. 


1. 0 < a < 1 at criticality 


For the critical model, we have A = 0. We study 
the two-point functions for the long-range hopping and 
pairing terms separately. 

Hopping .—In this case, the dispersion relation reads 

approximately Wa(g) = A\^\ -f « W A +\\q?~‘^'^)- 
Using Eq. (34), one can show that the correlation func¬ 
tion behaves asymptotically as 


GAR) 


/ llp>3-2a \ 


(40) 


T)cr(t, R) OC 


1 

2 



dq 

^a{q) 


sin(a;cr(g)t)e*'^''^ 


( RAq) 

[vA-^q) 


VAA) \ 
-Ra{q)) ' 
(37) 


For the sake of comparison, we first quote the two-point 
functions in the absence of long-range terms at or away 
from criticality. Within the short-range model, these 
functions are directly related to their counterparts in the 
case of the scalar field via the first terms in Eqs. (34) and 
(36). The correlation function at criticality, A = 0, is 


GsAR) - f GsAR) 


-X/R? 1/i? 
-1/i? XjR? 


(38) 

where, in the last step, we used GsR(i?) = GsR(i?) 
log R. Away from criticality, the correlation function falls 
off exponentially beyond the correlation length. Simi¬ 
larly, the response function at criticality is given by 


Note that we have not kept track of the precise coef¬ 
ficients. At long distances, the dominant power-law is 
given by 1/7?^“'’’. The corresponding exponent, 2 — u, 
is not consistent with the scaling dimension of the Fermi 
field, ['k] = 1/2; however, we have dropped a leading- 
order delta function 5(7?) which has the correct scaling 
dimension, but nevertheless is purely local. We remark 
that, even if one reads off the scaling dimension of the 
Fermi field from the exponent in Eq. (40), interaction 
terms still would be irrelevant, and our field theory up to 
the quadratic order in the fermionic field is exact in the 
RG sense. 

To find the response function, we shall focus on the 
region well outside the linear light cone t = R, where one 
can truncate the dispersion relation at the leading order 
as uJa{q) « Using Eq. (36), one finds 


0,^(7, 7?) 


idt Or 


-dR 

idt) 

b 

1 

1 0 


(^ 

0 -1 

7?" 

\r'^ 


7?- 


Dsr(?, R) 


idt 

-Or 


Or 

idt 


0(7 - 7?) - 5(7 - 7?), 


(39) 


The first term in this equation is obtained by acting with 
linear derivatives on the response function computed in 
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Eq. (16) with h„(s) = g 2 a{s) [the function is defined in 
Eq. (16)], while the second term is due to in Eq. (36), 
with ha defined as 

pOO 

hais) = / dq sin {q'^s) COS q . (41) 

Jo 


The derivatives in the above expression for D^- can be 
organized to cast it in terms of two separate scaling func¬ 
tions 


Ba{t,R) 


: ^ii( R'’) 


1 

i? 2 -a 


n2[ 





(42) 


where 


/iii(s) =ih'^{s) + ha{s) , 

hi 2 {s) = -(1 - (j)ha{s) - (Tsh'J^s) , (43) 

and the prime on ha indicates the first derivative. Thus 
all the terms in the response function feature non-linear 
causal regions with the relative scaling of time and space 
coordinates satisfying 


t ~ i?'", (44) 

consistent with the dynamic critical exponent z = a. 
The dominant contribution at long times and distances 
is given by the diagonal terms of Eq. (42) and is of the 
form 


with 


T)a{t,R) 



H,(s) 


/ hii{s) 0 \ 

\ 0 /l22(s)/ ’ 


(45) 


(46) 


Finally, at long distances and at a fixed time, one can 
see that the dominant term is {\/R)h{t/R"^) ~ 

[since h{s) ~ s for small s], again consistent with 
the Hastings-Koma bound [22] for the lattice model of 
fermions with long-range hopping. 

Pairing .—In this case, the dispersion reads uJa{q) = 
V|g|2- + |q|i+-« |g|-(l+i|g|i—). Using Eq. (34), the 
correlation function becomes 


Ga{R) 




(47) 


At long distances, the dominant power-law is given by 
l/R^~’^. We refer to the paragraph following Eq. (40) 
regarding the associated exponent in comparison with 
the scaling dimension of the Fermi field. 

Similar to the hopping case, in order to find the re¬ 
sponse function in the region well outside the linear light 
cone t = R, we can approximate the dispersion relation 


as u!a(q) « One can go through the same argument 
as above to obtain the asymptotic (at large distances and 
times) response function as 


T>a{t,R) 


1 

R 



where 


kii{s) = ih'J^s ), 
ki2{s) = ka{s), 


(48) 


(49) 


and ka is defined as 

poo 

ka{s) = / dgsin (g°’s) cos ((7-I-7r(T/2) . (50) 

Jo 

Again, we find a nonlinear causal behavior described by 
Eqs. (44) and (45) with the identification 


H,(s) 


/ fcii(s) A:i2(s)\ 
\-ki2{s) kii{s)J 


(51) 


Furthermore, the scaling with R at fixed t is consistent 
with the Hastings-Koma bound applied directly to the 
lattice model of fermions with long-range pairing (as can 
be seen by noting that, similar to ha{s), we have ka{s) ^ 
s for small s). 


2. a > 1 at criticality 

For a > 1, the dispersion relation can be approximated 
at small [gl as oja{q) ~ j^j- The correlation function is 
then given by 

r (R) f^( iq + Va{iq) 

J q \-iq + Va{-iq) -q^-Ra{q) 

In the case of long-range hopping, we find 

with a = min(2,cr). The dominant power-law at long 
distances is given by 1/R. In the case of long-range pair¬ 
ing, we find the same asymptotic dependence on R as in 
the short-range model, Eq. (38). 

The response function can be computed with the aid 
of Eq. (36), the first term of which simply consists of 
derivatives acting on Da{t, R) from the scalar case. The 
dispersion relation to the subleading order takes the form 
uja{q) « l^l -I- Ba\q\°'~^ with a = min(2CT, 2 -|- a) and 
some constant Ba- Therefore, the function Da{t,R) can 
be taken from Eq. (19) with the identification cr —>■ a = 
min(2cr, 2 -|- cr), that is Dait, R) = Da{t, R). Notice that 
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a > 2 for cr > 1, hence, D{t,R) in the first term of 
Eq. (36) exhibits a linear causal behavior when cr > 1; 
derivatives acting on this term should not change the 
causal behavior. Also, the second term in Eq. (36), 
Do- (t, R) , can be computed asymptotically at long times 
as (assuming ujcr{q) ~ |(z| and R > t) 

'Dcr{t,R)= j dq q'^~^ sm{qt) cos{qR) 

1 1 

"" {R-ty ~ {R + ty' 

Near the light cone t = R, we have Do.(t, i?) ~ 1/{R — 
ty at long times, indicating a linear causal behavior. 
Therefore, the full response function T>„ indeed gives rise 
to a linear light cone. Furthermore, at long distances 
and fixed time, Dcr(t,i?) is again consistent 

with the Hastings-Koma bound [22] for the corresponding 
lattice model. 

For long-range pairing, the dispersion relation to the 
subleading order in jgj is given by uja-{q) = jgj + Byqy~^ 
with a = min(2cr, 1 -|- cr) > 2 and some constant By 
once again a > 2 for cr > 1. One can argue, in a similar 
manner to the case of long-range hopping above, that the 
causal region is bounded by a linear light cone. 


3. All (7 > 0 away from criticality 


For long-range hopping, and setting A = 1 for con¬ 
venience, the dispersion relation is given by uia{q) = 
yi-\- q^ + jqj'^, which yields, via Eq. (34), 


GAR) 




(53) 


At long distances, the dominant power law is given by 
the diagonal contributions ^ 

For long-range pairing, the dispersion relation is given 
by UJa-{q) = y/l + q^ + \qy with a = min( 2 cr, 1 -|- cr); we 
find 


G^R) 




(54) 


in agreement with Ref. [40]. At long distances, the dom¬ 
inant power law comes from the off-diagonal contribu¬ 
tions ^ We note that, in both cases above, the 

correlation functions are described by the same asymp¬ 
totic power-law as that of the long-range couplings in the 
lattice Hamiltonian. While this appears to be a fairly 
generic feature of long-range interacting models, there 
are exceptions [41]. 

The response function is given by Eq. (36) and can be 
cast as 


R) ~ ±Dcr{t, R) + higher derivatives, (55) 


where those terms that include higher derivatives or 
higher factors of momenta in the integrand (D^) are ne¬ 
glected. Therefore, information about causality in the 
non-critical fermionic system is captured by D„ {t, R ), 
and is essentially identical to that of a scalar-field model 
with the same dispersion. Consequently, the conclu¬ 
sions of Sec. Ill A 3 should be immediately applicable 
here. Specifically, for long-range hopping, ui„{q) = 
•^1 + q'^ + IqI'^, and we find a non-linear light-cone for 
a < 1. On the other hand, for long-range pairing, 
uja{q) = \/l + q'^ + \q\°‘ with a = mm{2a, 1 -|- cr). The 
nonlinear light cone then arises for a < 1 , or equivalently 
a < 1 / 2 . 


V. Conclusion and Outlook 


In this work, we have studied the critical and near- 
critical properties of the long-range interacting TFIM 
and of a fermionic model with long-range hopping and 
pairing. We have identified the critical exponents char¬ 
acterizing the dynamic behavior as well as the two-point 
correlation functions of the system. For the TFIM, a 
nontrivial calculation gives the value of these exponents 
to the two-loop order. Field-theoretical arguments indi¬ 
cate that the fermionic model is exact (in the RG sense) 
at the quadratic level, from which we are able to cal¬ 
culate various critical exponents. At the critical point 
, we have argued that the causal behavior is identified 
by the dynamic critical exponent, and the response func¬ 
tions find a general scaling form, explicit and universal 
properties of which have been obtained. For both critical 
and noncritical models, we have derived the regimes of 
linear/non-linear light cones, mostly with the aid of the 
dynamic critical exponent. 

Our definition of the causal behavior and causal region 
is inspired by field theory; specifically, the linear causal 
region at criticality is identified with z = 1 , indicating 
relativistic dynamics. We have argued that values of 
z < 1 correspond to a sublinear causal region. It is shown 
that the causal region defined by the dynamic critical 
exponent indeed characterizes the local maxima of the 
response function, suggesting a close connection to the 
speed of information propagation. We have also shown 
that the response function can generally be expressed in 
terms of scaling functions up to multiplicative decaying 
power laws. Recalling that, in a local quench, the sig¬ 
nal typically decays with distance, we find our definition 
of the causal region based on the scaling function quite 
natural. However, we stress that this definition need not 
produce the same light cone that would be obtained by 
setting the tightest possible Lieb-Robinson-type bound 
equal to a constant, a subject that merits further inves¬ 
tigation. If the two agree at criticality, where our cal¬ 
culations predict the strongest deviations from a linear 
light cone, we conjecture that the TFIM and the Fermion 
Model exhibit a linear causal behavior for I -|- cr ^ 3 
and 1 -f cr > 2, respectively. This conjecture is particu- 
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larly valuable in light of the fact that even the tightest 
available Lieb-Robinson-type bound for long-range inter¬ 
acting systems cannot exclude the existence of a nonlin¬ 
ear light cone at any finite tr, no matter how large [28]. 
It would also be worthwhile to explore possible connec¬ 
tions between this conjecture and the proof that the light 
cone cannot be logarithmic, and is at worst algebraic, for 
1 -bcr > 2 [28], 

The powerful tools of field theory used in this paper 
can be employed to study a range of interesting, and 
experimentally relevant, long-range interacting systems 
[16, 53-56], such as a huge variety of spin-1/2 [57-60], 
spin-1 [41, 61, 62], and higher-spin [63, 64] models, gen¬ 
eralized Hubbard [63, 65, 66] and t-J models [58, 59], 
and spin-boson problems [67], among many others, in 
one or more spatial dimensions. In general, these mod¬ 
els exhibit new universal behavior not captured by stan¬ 
dard long-range interacting classical models, since the 
quantum-to-classical mapping generates classical models 
with long-range interactions in all spatial directions ex¬ 
cept the one corresponding to the imaginary time dimen¬ 
sion of the quantum model [39] . Another interesting, but 
more complicated, direction is to study the consequences 
of a global quench as opposed to a local perturbation. 
However, in this case, it is not at all obvious that field- 
theorertic considerations will be valid given the extensive 
energy imparted by the quench (though see Ref. [68] for 
related studies of weak global quenches). Finally, this en¬ 
tire line of speculation is intimately related to the growth 
and propagation of entanglement in systems with long- 
range interactions [22, 23, 69-71], a subject for which our 
understanding is far from complete. 
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A. Two-point functions in quadratic scalar and 
fermionic field theories 


Correlation function. —The time-ordered correla¬ 
tion function in real time, G{t,R), is related to Q (in 
imaginary time and real space) as [72] 

iG(t, R) = Q{it + sgn{t)0^, R). (A2) 

Thus, 

iG{t ^ 0, R) = G{0, R) (x j dujdqG{uj,q)e^‘^^, 

which leads to Eq. (12) for the correlation function at 
equal times. 

Response function. —The causal response function 
D{t, R) is related to G via an analytic continuation in 
frequency space as iuj ^ uj + iO''" [72] , that is 

D{u;,q) = = (a, + i0+)2 -w(g)2 ' 

(A3) 

This equation, cast in real space and time, leads to 
Eq. (13). 


2. Fermionic field 


For the fermionic field, the imaginary-time-ordered 
correlation function is obtained by inverting the 2x2 
matrix in Eq. (33), 


GF{u},q) 


fiuj + A - q'^ 
\ -iq 



0 (w,(?) 


-b 


( nq) 

\V{-iq) 


V{iq) \ 

-n{q)) 


0(w,g), 


(A4) 


with G from Eq. (Al) (with uj{q) specific to the fermionic 
theory and the subscript a dropped). The first line of this 
equation can be cast in terms of the function G and its 
spatial/tempoal derivatives, while the second line can¬ 
not. A similar analytical continuation described for the 
scalar case leads to Eqs. (34) and (36), wherein the tilde 
functions correspond to the second line of Eq. (A4). 


In this appendix, we derive the analytical expressions 

for the two-point functions in a quadratic field theory in ^ rq for the cj)* model with long-range 

terms of the dispersion relation U){q). interactions in the spatial direction 


1. Scalar field 


For a quadratic scalar field theory, the imaginary-time- 
ordered correlation function in Fourier space is given by 

9) = 2 2 ( \2 ’ 

a;2 + u]{qy 

with ijj{q) the dispersion relation. 


In this appendix, we perform the RG calculation in¬ 
cluding the two-loop diagrams for the (j)'^ model in 1-1-1 
dimensions with long-range interactions along the spatial 
direction. Although we are interested in d = 1 spatial di¬ 
mension, we generalize to a d-dimensional space in order 
to expand around the upper critical dimension d = d„. 
We first cast the action in imaginary frequency and mo¬ 
mentum space as 
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H = 


+ u 


r duj g + Aw'^ + B^q<^ 

(27r)^ J ^ 2 

r d‘^qid'^q2d'^q3 f dujidui2dui 

(27r)3d J {2tt)^ 


</'(w,q)(/>(-a;,-q) 

qi)(/i(w2, q 2 )<?i(a; 3 , q 3 ) 0 (-a;i - W3 - ^3, -qi - q2 - qs), 


(Bl) 


where q = |q|; note that we are using a different nor¬ 
malization (a factor of 1 / 2 ) for the quadratic terms in (j) 
for convenience. (In a slight abuse of notation, we use 
the same coefficients A, B^, and g defined in the action 
cast in real space and time coordinates.) In order to 
compute loop diagrams, we rely upon the momentum- 
shell RG [47]. However, similar to Ref. [73], we do not 
find it necessary to introduce a cutoff in frequencies, and 
will integrate over oj S (— 00 , 00 ), which is free of di¬ 
vergences. Without long-range interactions, time and 
space coordinates appear on equal footing in the action 
+ u^(V^)^ for some velocity u], and long wave¬ 
length physics becomes ‘Lorentz’ invariant, which yields 
z = 1 exactly. For long-range interactions, on the other 
hand, there is no such symmetry, and it is natural to 
treat space and time coordinates differently as prescribed 
in the main text. 

We first compute the one-loop correction (depicted in 
Fig. 5) to the parameter g in the action; we find 


'4\ r d'^q r°° duj 1 

2/ ^ JA/b<\q\<A (27r)^ y_oo 27r £< + Auj"^ + B^q'^ 
'4\ u r d^q 1 

2/ 2^/A JA/b<\q\<A (27r)‘^ Ve + B^q<^ ’ 


where we have used the identity ^ 

The constant A is the momentum cutoff in the model, and 
the integral is over the momentum shell A/b < |q| < A. 
We specialize to an infinitesimal rescaling b = 1 + 61 and 
expand everything to first order in 61. The RG flow of g, 
combined with its scaling dimension in Eq. (26), is then 
given by 


dg 


{z + d — 2a)g + 


6KdA‘^u 
\/ A-^/g BcjA^ 


0{u^), (B2) 


where we have used the identity “ 

— A— The RG flow for u, combined with Eq. (26), 
iVAB^ 

reads 


du , , ■ X 

— = (z-|-a —4a)M — 
dl 


9KdAV 


VA{g + B^A’^)3/2 


-0{u^). (B3) 


At the one-loop order, neither A nor B^r are renormalized, 
and thus their RG flow is simply given by the engineered 
scaling dimensions in Eq. (26), 


—^ = {z - a + d - 2a)B^ + O(u ^), 
dl 
d A 

\ = [-z + d-2a)A + 0{u'^). (B4) 

dl 


To this order, the above set of equations determine a = 
{d — a 12) 12 and z = a 12 (identical to the mean-field 
values). The RG equations (B2) and (B3) for g and u 
then take the form 


dg QKdA'^u 

dl ^ VAy'g + B„A<^ ’ 

du 9KdA'^u^ 

-— = eu -;=-— . (B5) 

dl VA(p + R,A-)3/2 

The fixed point is obtained by dg/dl = du/dl = 0, and 
yields the fixed-point values u = and g = g* given by 


— 


2B^A'^ 

3tT 


e, 


u* 


y/ABl'^ 


(B 6 ) 


where we have kept the dependence on e only to the first 
order, thereby replacing other appearances of the dimen¬ 
sion d hy du- Now the RG flow near the critical point 
can be linearized as 


where Kd = Sd/{2Tr)'^, with Sd the surface area of a 
sphere in d dimensions. 

Next we compute the one-loop RG flow of the interac¬ 
tion vertex. In computing the corresponding diagram, we 
drop the dependence on external momenta and frequen¬ 
cies which otherwise produce derivatives in the nonlinear 
terms, and therefore can be neglected [50]. We obtain, 
for the one-loop diagram, 

/ 4 \ r ^ r- d^ _ 1 _ 

\2/ V2/ JA/b<\q\<A {2n)‘^ 27r (^ -|- -|- Baq"^)^ 

_ r d^q 1 

VA yA/6<|q|<A (27r)‘^ _|_ i3^g'T)3/2 


^ (6g\ ^ (a-e/3 • • • W 6g\ 

dl \6uJ y fA(e^) —ej \6uJ ’ 


(B7) 


where 6g and 6u are deviations from their correspond¬ 
ing fixed-point values. The top element of the second 
column is not computed as it is not necessary for com¬ 
puting eigenvalues [50]. The above eigenvalue equation 
can be used to compute various critical exponents. For 
the critical exponent defined via the divergence of the 
correlation length near the critical point, ^ ^ {6g)~^, we 
find 


1 

V 


= a — 


e 

3 


+ 0{e^), 


(B 8 ) 
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the critical exponent characterizing the divergence of sus¬ 
ceptibility, X ~ {Sg)~^, is given by 

7=1 + ^+ O(e^), (B9) 

and, the critical exponent rj characterizing the anomalous 
decay of the correlation function, G{R) ~ be¬ 

comes 

r]=l-a/2 + 0{e'^), (BIO) 

to the first order in e-expansion. Together with the dy¬ 
namic exponent z, the above exponents satisfy the gen¬ 


eral identity 


j = v{z - ri + 1) . (Bll) 


Next we go to the two-loop order to find corrections to 
the exponents z and rj. We first note that knowledge of 
the fixed-point values of g and u to first order in e suffices 
to find the above exponents at the order. To this end, 
we need to compute the two-loop diagram in Fig. 5; we 
find 


2 f d‘^qd‘^p r°° dJduj" 

[yvr A/.<iqi,ipi,iwi<A7^y-oo7^ 

_ y^ _ 

g +Auj''^ + g + + B^p'^ g + A{oj'+ to” - uj^ + Bcr\q +p - k]’^ ' 

where w = q -|- p — k, and k and uj are the external momentum and frequency, respectively. This expression possibly 
gives corrections to the frequency- and momentum-dependent terms, i.e. the coefficients A and in the quadratic 
part of the action (Bl). To find such corrections, we must expand the above integral for small frequencies and 
momenta keeping in mind that |a;| ^ |a;'|, \l>j"\ and |k| ^ |q|, |p|. A Taylor expansion in the former will produce 
only analytical terms, and thus terms such as q'^ cannot be generated in this procedure. This implies that Bn is not 
renormalized, and leads to the non-renormalization condition 

z — a + d— 2a = {). (B13) 

This is an exact equation, and is not affected by higher orders of perturbation theory. Combined with a = {d— l+r])/2, 
this equation yields the exact relation between the critical exponents rj and z in Eq. (28). In light of this identity, we 
find that Eq. (Bll) takes a particularly simple form, 


7 = (71/ . 


(B14) 


Next, we consider the renormalization of the term in the action. We should expand the expression (B12) for the 
two-loop diagram to second order in w, 


OC 


, 2 , ,0 


, 48u^ 2 

)- 


d‘^"qd‘^“p 


1 


1 


4^3 7l/6<|q|,|p|.|w|<l (277)2^^- q^Gp'^Gw’^G 


3 ’ 


(B15) 


where we have assumed that g is small, justified by the fact that g ^ e at the fixed point. The first term, being 
constant in w, renormalizes g to order O(e^), but is inconsequential in determining the renormalization of A at this 
order. With the knowledge that u* ^ e, we have also replaced d by the upper critical dimension du = 3(t/ 2 in order 
to keep only the leading order in e. In computing the second term in Eq. (B15), we have used the identity 


r°° dJduj" 1 1 1 

J_^ (27r)2 Auj''^ + a? -I- iP Ay + uj" — ojY + 

a + b + c y , 4 n 

AAabc [(a -I- 6 -b c)^ -|- Acu^] 4a6c(a + b + y 

Also, in the momentum integral, we have rescaled all momenta by A; powers of A from the integral measure and 
the integrand cancel each other out. It is worth pointing out that, unlike the one-loop integral, one cannot consider 
an infinitesimal rescaling. In fact, it is convenient to consider a finite b, and reiterate the renormalization group 
transformation n times. To evaluate the integral in Eq. (B15), we note that a rigid rotation of all three momenta 
q, p, and w does not change the integral. We thus fix, say, q = (^z in the direction z, and cast the integral over 
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q as / j 2 ^ ~ I dqq'^'^ The integral over p can be then cast in polar coordinates, in du dimensions, as 

I { 2 TT)du = /o^/ dpp‘^'^~^. In these coordinates, the integration limits are given by 1/b < q,p < 1, 

and 


1/6^ - q^ -p^ 
2qp 


< cos 9 < 


I - q^ -p^ 
2qp 


A change of variable p = qy allows us to write the momentum integral in Eq. (B15) as (with du = 3(t/2) 


^3(T/2-f^3cr/2-l 


r dq 
'l/b 9 , 




dy / dO 


sin0) 


■ia/2-2 


(1 + 2/2 + 2j/ cos Oy/'^ [l + 2/‘^/2 _|_ (1 -|_ y2 _|_ 2 y cos 0)°’/^]' 


(B16) 


where the limits of integration take the form l/{bq) < y < 1/q and a more complicated expression for cos0. However, 
since we are only interested in the leading logarithmic dependence, log b, we can simply replace the limits by 0 < // < oo 
and 0 < 9 < TT without encountering any divergences. For example, for // —t 0, the integrand goes as y'^~^, while, for 
2 / —)■ 00 , it goes as 1/y'^'^^, thereby being perfectly convergent in both cases for cr > 0. This simplification allows one 
to directly compute the integral over q to obtain a factor of log h. The remaining integral over y and 9 then directly 
contributes to the two-loop diagram as 

const — (log 6) o(tT)e^ , (B17) 


where 


o(cr) = 


4r(3cr/2) 


27y^r(3cr/4- 1/2) 


dy / d9 


2/'"-i(sin6l)3<"/2- 


0 Jo (1-I- 2 / 2 -f 2//cos [l-I- 2 /'^^^ + (1 + + 22 /cos 0)°’/"^]' 


(B18) 


Note that in computing this expression we have used the fixed-point value of m* given by Eq. (B6). We have also 
verified numerically that the exact evaluation of the integral in Eq. (B16) almost exactly reproduces the log 6 term in 
Eq. (B17). The recursion relation for A„ is then given by 

A„+i = A„ 6'^-2- [1 + 2(log b) ^(a)e2] « , (B19) 

where the exponent of b in the first equality is obtained from the scaling dimension in Eq. (26) [or, equivalently, 
Eq. (B4)] combined with the non-renormalization condition, Eq. (B13), while the expression in the bracket gives the 
two-loop correction. In the last equality, we have exponentiated the e-dependent term justified by our perturbative 
treatment in e. The above recursion relation goes to a finite fixed point A* if the exponent z = a/2 + Az is given by 
Eq. (29). Similarly, the exponent f] is computed via the exact relation (28). For the special case of cr = 1, the above 
exponents have been computed in Ref. [49], 


1 (N + 2) (l2-7r2) „ 

2 16(iV-b8)2 


Oiey, 


for an iV-component Landau-Ginzburg field theory in d < 2 spatial dimensions, where e = 2 —d. With d = 1 dimension 
(e = 1) and = 1 component, we find that the above result is consistent with Eq. (29) where we should substitute 
a = I and e = 3a/2 — 1 = 1/2. 

To find the behavior of our model away from the critical point, we note that the renormalization group does not 
change the power laws obtained at the quadratic order, specifically, logarithms such as the one in Eq. (B17) do not 
appear. 

Finally, we make the comparison between our results up to the one-loop order with Ref. [39]. A different renor¬ 
malization scheme has been adopted by the authors of Ref. [39] where they have also integrated over a frequency 
shell, while we have integrated over the whole range of w at each step of RG. Therefore, their RG equations seem 
different from ours in Eq. (B5); however, we have checked that the linearized equation near the fixed point, Eq. (B7), 
is indeed the same in both cases. While our exponent iz is different from the one in Ref. [39], we believe that the 
reason is a typographical error in Ref. [39] . Furthermore, the exponent rj in Ref. [39] associated with the spatial cor¬ 
relation function seems—although not explicitly stated—to be different from our definition: the definition of Ref. [39] 
is G{R) ^ , compared to our definition G{R) ~ 1/7?'^“^+''. 
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